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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
eoRPORAT ibe NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONS . 


IBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH 


L 
S NOT SUPPLIED BY DIGITAL. 
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of ACaLity Language support Library - user callable 
ABSTRACT: 
Integer word base to integer word power. 


Integer overflow can occur if the result exceeds a word. 
Undefined exponentation can occur if base is 0 and power is 0 or negative. 
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VERSION: 0 
HISTORY: 
AUTHOR: 
Thomas N. Hastings, 5-May-77: Version 0 


MODIFIED BY: SUSAN HUBBARD AZIBERT 
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OTS$POW!1 1 = INTEGE ER*2 ** INTEGER*2 power routine 16-SEP-19 $ 06: iT ¢] AX/VMS Macro V04-00 
a 806 ISTORY ; Detailed Current Edit History 6=-SEP-1984 7 CMTHRTL.SRCJOTSPOW!1.MAR; 1 
-SBTTL HISTORY ; Detailed Current Edit History 
4 ; Edit History for Version 0 of OTS$POWI! 
5 ; version 04 = changed module name to forpow!! 
§ i; version ; = changed error handler from ATHSERROR to MTHS$SERROR 
3 arg ree 9 - changee error handler to MTH$$S 
0 3 3 O-1 fixed bug in case instruction. JMT wheb-78 
0 ; 1-001 Update version number and copyr ht notices goat Pe aati 
000 60 ; 1~308 = Change | x UNDEXP to MTHS$K_UNDEXP. 
000 61 ; 1-003 = Add o the PSECT directive. JBS BS eco 
000 6¢ : 1-004 - Use 32-bit addresses to reter te te. gternals. JBS 28-JAN-1979 
000 63 ; 1-005 = Declare externals. SBL 17-May-1 
0000 64 ; 1-006 = Use general mode addressing. "Sar 30-Nov-1981 
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-SBTTL DECLARATIONS 


INCLUDE FILES: 


EXTERNAL SYMBOLS: 


-DSABL GBL 
-EXTRN MTHSK_UNDEXP 
eEXTRN MTHSSSTGNAL ; Math error routine 
80 3; MACROS: 
00 3 
00 
00 ; 
38 ; EQUATED SYMBOLS: 
00000004 000 base = 4 ; base input formal - by-value 
00000008 33 exp = 8 3 exponent intpu formal - by-value 


; OWN STORAGE: 


PSECT DECLARATIONS: 
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-PSECT _OTSSCODE PIC,SHR,LONG, EXE, NOWRT 
3; program section for OTS$ code 
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-SBTTL OTSSPOWII - Word to power word giving word result 


p++ 
; FUNCTIONAL DESCRIPTION: 


Signed word result = signed word base ** signed word exponent 
The signed word result 


exponent 


> 


AAAAAA 


IMPLICIT INPUTS: 
NONE 


FUNCTION VALUE: 
Word integer base ** exponent 

SIDE EFFECTS: 
SIGNALS SSS$_ARITH with integer overflow hardware code if 


Ow. 
NDEXP (82 = ' UNDEFINED EXPONENTATION') if 
base is 0 and exponent is 0 or negative. 


0 


Ooooooeo Goo 


CALLING SEQUENCE: 
Power .ww.v 

INPUT PARAMETERS: 
NONE 


OUTPUT PARAMETERS: 
NONE 


IMPLICIT OUTPUTS: 
NONE 


nteger overfl 
SIGNALS MTH$_U 


and even 
and odd 


Page 


S given by: 

result 

product (base * 2**i) where i is each 
non-zero bit position in exponent 

1 

greet rmes exponentation 


0 

1 

Undefined exponentation 
-1 

1 


Integer overflow can occur. 
Undefined exponentiation occurs ‘f base is 0 and 
exponent is 


or negative. 


OTS$POWI! (base.rw.v, 


exponent.rw.v) 
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1986 99:98:87 


ENTRY OTSSPOWII, “M<IV, R2> 
MOVL 


H 
INTEGER*2 power rou 
rd to power word giv 
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: Exponent > 0. 
; Scan each exponent bit from ri 
: For each l-bit in exponent, mu 
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; R1 = ba 
; undefin 
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exponent =< 0. > 
Result is given by the following table: 
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; enable integer overflow 
RO = i 


initia 


= exponent 
xponent =< 0 


; R 
; branch if e 


R1 = base 


result 


ht, squaring base each time thru loo 
tiply current base into partial resu 


branch if Supqgens is odd 


R2 = exponent/ 
R1 = current 


ower of base 


integer overflow will trap 
ARIT 


and SIGNAL SS 
loop if next bit in exponent is 0 


next bit in exponent is a 
RO = new partial result 


integer ove 
and SIGNAL 


R2 = exponent/ 
f more exponent bits are 1 
> = base ** exp 


pep 
re rn 
RO<31:16> = 


0 too, else return 1 


RO<15:0 


se 

ed if base 
; return with result 
; since base ** 0 


iLL tra 
RITH 


; branch if exponent < 0 


= 0 too 


or -1 depending on exponent being even or odd 
pndet ines exponentation 
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1856 OTSSPOWI1 = Word to power word stvine wo pease 7 99:98:3) EATMRTL. SREJOTSPOW! 1 .MAR: 1 . (4) | 
15 EXPLSS: 
O02 -FFFF 8F 04 AC _ AF 18 ASEW base(AP), #-1, #2 3 fas on value of base 
0009" 17 108: .WORD MINUS1-16$ ; C-1J: return RO = -1 or 1 depending 
) 1 i; ON_exponent being odd or even 
OO ape 8 1 WORD oo | in | ha ; C0]: Undefined exponentation 
OF' 0030 920 {WORD POWIIX-T08 ; (+1); return RO = 1 | 
50 BS F 1 CLAL R 3; [< -1 or > #1): return RO = 0 
4 ry § RET 
7 4 MINUS1: 
03 52 «£9 4 5 LBC R2, POWILIX ; if exponent is even, return RO = 1 
50 01 fF Be § MNEGL #1, RO ; else return RO = -1 
4 $8 POWIIX: RET 3 return 
04 8 
049 9 3+ 
Bpee ? 3; Undefined exponentation error = 0**0 or 0**(-n) 
049 $ 
049 UNDEFINED: 
50 D& 0049 4 LRL RO 3 return result = 0 if error 
7E 00° 8F 9A 0048 5 MOVZBL #MTHSK_UNDEXP, -(SP) ; FORTRAN error # 
00000000'GF 01 FB O04F 4 CALLS #1, G*ATHSSSIGNAL 3 convert to 32-bit condition code 
0056 7 ; and SIGNAL FORS_UNDEXP 
04 9038 8 RET 
005 39 
0057 40 . END 


§ 3 
= INTEGER*2 ** INTEGER*2 power routine 


OTSSPOWI1I 16- 4 AX/VMS 0 V04-00 o 
Symbol table greeks 99:98:8) FANANES Bacto yoRe OO ans '29e 
BASE 3H 
EXP 
EXPLEQ 1 
Exel es ahs 1 
ad et 1 
MTHSSSIGNAL eereeres x 0 
MTHSK_UNDEXP 90000000 x 0 
s$POw1 1 £008 00 RG 1 
PARTIAL £000 iD R 1 
I 48 R 1 
SQUAR 00 & 1 
QUAR1 0 7 R 1 
UNDEF INED 49 R 1 
eee ee Ree me ee + 
! Psect synopsis : 
PSECT name Allocation PSECT No. Attributes 
ABS 00000000 <( 0.) 00 ¢ O.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
OTSSCODE 00000057 <¢ 87.) O01 ¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
powwow ewer ree ene nnn en cmrone 
! ; Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 1 00:00: 8.06 00:00:01.01 
zonnend processing 128 00:00: 8°32 88:60:03" H 
Symbol table sort 50 00:0 : 0:61 3:81 8 
Symbol table output : 00:00: 0:03 $5 :00:0 -08 
Psect synopsis output 0:00:00.0 00:00:00.1 
Cross-reference output § B92 82: 9-99 Bp #28 +88 -08 
Assembler run totals 28 0:00:01.76 00:00:09.60 


e working set Limit was 900 pages. 
35%9 bytes {8 pages) of virtual memory were used to buffer th 
There were pages of symbol table space St iecetee to hold 13 non-local and 1 Local symbols. 
40 source Lines were read in Pass 1, produci t. 1 object records in Pass 2. 
pages of virtual memory were used to define 0 macros. 


intermediate code. 


ee eee see ener newer ere eer eceee $ 


Macro library name 
.$255$DUA28: (SYSLIBISTARLET.MLB; 2 0 
0 GETS were required to define 0 macros. 


There were no errors, warnings or information messages. 
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MACRO/ENABLE=SUPPRESSION/D1SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOW11/OBJ=OBJ$:O0TSPOW11 MSRC$:OTSPOW! 1/UPDATE=(ENH$:O0TSPOW! 1) 
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